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The interplay between the ground state energy of the generalized Bernasconi model to multi-phase, 
and the minimal value of the maximal autocorrelation function, Cmax = raaxK \Ck\, K — 1, ..N — 1, 
is examined analytically and the main results are: (a) The minimal value of miuiv Cmax is 0.435%/^ 
significantly smaller than the typical value for random sequences 0{^/\og Ny/N). (b) minjv Cm a ^ 
over all sequences of length A'' is obtained in an energy which is about 30% above the ground- 
state energy of the generalized Bernasconi model, independent of the number of phases m. (c) 
The maximal merit factor Fmax grows linearly with m. (d) For a given A'^, miUjV Cmax ~ \/ N /ttl 
indicating that for m = N, miujv Cmax = 1, i-c. a Barker code exits. The analytical results are 
confirmed by simulations. 



In many applications of communication science 
as well as a variety of other fields, it is necessary to find 
sequences with low autocorrelation functions. Some of 
these applications utilize the pulse compression feature 
of the low autocorrelated sequences to obtain high res- 
olution in radars and sonars. In other applications, the 
shifts of such periodic sequences can be used to identify 
users in multi-user systems. Due to their importance, low 
autocorrelated sequences have evoked a wide spread in- 
terest accompanied by the development of various meth- 
ods for constructing such sequences 

In order to construct a binary sequence S = (si, sjv) 
with low off peak autocorrelations, one has to define 
which quantity has to be minimized. The different appli- 
cations divide the low autocorrelated sequences into two 
types in correspondence with the quantity which they 
minimize. The first kind of low autocorrelated sequences 
minimizes the Hamiltonian of the Bernasconi model |^|^] 
which is given by: 
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where for the case of non-periodic boundary conditions, 
which is at the center of our study. 
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Note that for random 



and 

Ck = = l ^j^((j+K-l) mod Af) + 1 

sequences the average value of H in the non-periodic case 
is whereas for low autocorrelated sequences this 

energy is reduced by a merit factor F > 1 to N'^/2F. 
However, there are applications for which the maximal 
off peak autocorrelation, Cmax = uiax^f |Cif|, has to 
be minimized. The second kind of low autocorrelated 



sequences are the solutions of this minimization prob- 
lem. Note that for random is typically 



0{^/\ogN V N) . For a sequence of length N, the maximal 
possible ratio between the peak, |Co| = N , and the maxi- 
mal off-peak autocorrelations , | Ck \ with = 1 , . . . , iV — 1 
is iV/1. The only known binary sequences with this ra- 
tio arc the Barker codes of length 2, 3, 4, 5, 7, 11 and 
13. Obviously, the Barker sequences, when they exist, 
furnish a minimum for the two minimization problems. 
However, Turyn |^] has shown that no other binary codes 
such as this exist for any length less than 144 or for odd 
length greater than 13. An exact solution for these two 
minimization problems is known only for systems (up to 
N = 59) which are small enough to permit an effective 
exhaustive search [^,|| . Extrapolation of the ground state 
energies which were found for small systems using an ex- 
haustive search indicates Fmax = limAr^oo ^ — 8.5 
Moreover, F^ax was conjectured by Golay |l^,|ll| to be 
bounded from above by 12.324. 

The following questions regarding low autocorrelated 
sequences are still open and are at the center of our study: 
(a) Do the sequences that minimize Cmax minimize the 
energy of the Bernasconi model as well? (b) In case that 
the two minimization problems are not equivalent, how 
far are the energy values of the sequences which minimize 
Cmax from the values of the ground state energy? (c) 
Exhaustive search of sequences with length up to iV = 59 
show that the degeneracy of the ground state is of 0(1) 
and is bounded from above by 8. Is the degeneracy of the 
states which minimize Cmax is of 0(1) as well? (d) In this 
context, which of the two minimal quantities is easier to 
approach? Additional interesting questions arise when 
multi-phase sequences, whose terms are complex m-th 
roots of 1 for m > 2, are considered. For such sequences 



si = exp(« — /) 
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where 1=1,..., m and the correlations are defined as 
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where Sj+x the conjugate of Sj+x- The information 
about muhi-phase low autocorrelated sequences is more 
hmited, mainly since the configuration space grows like 
which makes the exhaustive search ineffective even 
for very small sequences. Hence, not much is known 
about the influence of the number of phases both on the 
maximal value of F and on the minimal value of Cmax ■ 
Is it possible to attain better solutions by increasing the 
number of phases? 

In this study we suggest an analytical technique to 
calculate the quantities of interest, namely, max^r F and 
minff Cmax for both binary and multi-phase sequences. 
An obvious lower bound of minjv,_F Cmax over all se- 
quences of length N and a merit factor F is obtained 
by assuming that all \Ck\ equal N / {2F). Neverthe- 
less it is tempting to investigate the relations between F 
and minjv,_F Cmax and to measure more accurately the 
optimal value miuAr Cmax over all sequences of length 
N similarly to the aforementioned upper bound given 
for Fmax 10 • For our analytical study we consider the 
non-periodic Bernasconi model. We approximate the se- 
quences to be random and the autocorrelation functions 
to be independent variables with the following Gaussian 
distribution 
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Using Eq. (||), the probability PpiCmax) that the au- 
tocorrelations of a sequence with a merit factor F = 
N'^/2H have an upper bound Cmax such that \Ck\ < 
Cmax is given by 
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Since there are 2^ distinct binary sequences of length 
iV, it is necessary that Pp > 2~^ in order for a sequence 
with the corresponding features, namely a merit factor 
F and an upper bound Cmax, to exist. Hence, by equat- 
ing PpiCmax) to 2~^ one can find the minimal upper 
bound among the upper bounds of all se- 

quences with a merit factor F. Moreover, the Gaussian 
distributions of Ck imply that mhiN^p Cmax must be of 
0{\fN) in order for PpiCmax) to be equal 2~^. Assign- 
ing minjv^i? Cmax = B{F)^fN and inserting the integral 
representation of 5 function in Eq. (||) , the saddle point 
method can be used to obtain the following set of equa- 
tions 
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where g{y, A) 
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Solving numerically this 



2(l-a) 

set of equations, it turns out that a solution exists only 
for a bounded region of F, where .215 < F < 12.324. 
These two limit values of F are alternatively obtained 
by taking the limits of the integration in Eq. (^) to in- 
finity. The existence of a lower bound is a consequence 
of the fact that as the energy becomes smaller the al- 
lowed values of the correlations decrease and the proba- 
bility for such sequences (Eq. (^) decreases as well. As 
was pointed out above, when this probability drops be- 
low 2~^, there are no more sequences with the required 
energy. Similarly, for high energy values, high values of 
correlations are required and the corresponding proba- 
bilities become smaller again. The analysis of B{F) in- 
side the allowed region reveals that the minimal value 
of B{Fc) = 0.435 is obtained for F^ = 8.839 whereas the 
maximal merit factor F — 12.324 corresponds to a higher 
value of B{F) = 1. Thus, the merit factor of sequences 
which minimizes Cmax is around 71% lower then Fmax- 
Exhaustive search results in the region N S [15, 28] indi- 
cate that the average ratio between the two merit factors, 
<<< F{mmCmax) » /Fmax >N, IS ~ 0.6771 where 
<< >> denotes the average over all sequences which 
minimize Cmax whereas < > at denotes the average over 
the different sequences length N. The deviation from 
the analytical result ~ 0.717 is attributed to finite size 
effects. Fig. |l| shows the behavior of B{F) derived from 
the numerical solution of Eqs. (0) and (^). The behav- 
ior of Pp{Cmax) divides this graph into three regimes. 
Starting with the large F regime, there is only one se- 
quence whose Cmax cquals B{F)\/N . This feature holds 
until Fc = 8.839 for which the minimal B{F) = 0.435 is 
obtained. Further decrease of F results in an increasing 
number of sequences with an upper bound of B{F)^/N. 
This increment terminates at F — 0.8 and for smaller 
values of F the number of sequences reduces to 1 as 
approaches 0.215. In order to compare the analytical be- 
haviour of B to that obtained from simulations, it is nec- 
essary to find a way to circumvent the mismatch between 
the analytical values of F and those obtained from sim- 
ulations. The rescaling of F by Fmax yields a parameter 
F/Fmax whose range of values, [0, 1], does not depend on 
a specific realization. This allows a comparison between 
analytical results and simulations. Results of an exhaus- 
tive search over sequences of length = 32 show that 
the behavior of B as a function of F/Fmax resembles the 
analytical prediction (Fig. 0). 
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Moreover, the existence of the three aforementioned re- 
gions of PpiCmax) is confirmed by counting the number 
of sequences with the minimal upper bound B for each 
of the micro-canonical ensembles with a merit factor F 
and length N = 32. The histogram is demonstrated in 
the inset of Fig. ^, where the merit factor is rescaled 
to F/Fmax- In the small F regime, the average num- 
ber of sequences increases with F until F/Fmax ~ 0.2 
from which the number of sequences decreases up to the 
plateau which starts at F/F,nax ^ 0.83. 
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FIG. 2. B vs. F/Fmax- The solid line represents the ana- 
lytical results while the circles represent the exhaustive search 
results for A'^ = 32. Note that the analytical Fmax is 12.324 
whereas for A*' = 32 Fmax = 8. Inset: a histogram describing 
the logarithm of the number of sequences with the minimal 
upper bound B for each of the micro-canonical ensembles with 
a merit factor F as was derived from an exhaustive search for 
iV = 32. 

Now we turn to study the influence of the number of 
phases both on the maximal value of F and on the mini- 
mal value of Cmax ■ We assign Dk and Ek to be indepen- 
dent Gaussian variables which represent the real and the 
imaginary parts of Ck respectively. In case that m = 2^ 
where p is an integer number greater than 1, these two 
parts have the same probability distribution 
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Under the same assumptions that have been used in the 
derivation of the binary case, Eq. becomes 
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FIG. 1. B{F) vs. F. This graph is composed of three 
regions corresponding to the three different behaviours of 
PpiCmax)- In the dashed region 0.215 < F < 0.8, the frac- 
tion of sequences PpiCmax) is getting larger as F increases, 
while in the dotted-line region, 0.8 < F < F^, this fraction 
decreases to 2~^ as F approaches Fc- For F > Fc (solid-line) 
the fraction of sequences with the minimal value BiF)\/N is 
constant and equals 2~^. 
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Similarly to the binary case, if Pp''{Cmax) = m ^ then 
B{F)^/N is the minimal value of Cmax for which there 
is a sequence with a merit factor F. The same procedure 
that has been used to calculate B{F) for the binary case 
results in the following set of equations 



A 



1 



2F A 



((1- A) ln(l-A) + A) 



-2 / \ncrf{B{F)g{y,Xj)dy = ln{m). 
Jo 



(11) 



2B{F) 



1 1 

2F^X^ 
exp{~BiF)^g{y,X)^) 
TT Jo eviiBiF)giy,X))giy,X) 



(ln(l - A) + A) + 



dy = 0, 



(12) 



where g{y, A) 



l-A(l-y) 



(^i_y^ — ■ Solving numerically these 

two equations, the behaviors of minjv Cmaa; / and 
Fmax as a function of m are obtained. Fig. ^ shows 
that miuAT Cmaxj^TN is proportional to 1/ y/rn with slight 
deviations for small ni. This relation between BiFc) = 
min^r Crnax/^fN and m implies that such a relation holds 
for B{Fmax) as well. Since the variance of C\ equals 
N — K, their typical values drop linearly with K. Hence, 
H — N'^/2F is approximated by a sum of an algebric 
series of iV — 1 terms, C^, with an upper bound, C^ax^ 



3 



of 0{N/m). Assuming that C|- is homogeneously dis- 
tributed between to Cmaa; yields a linear increment of 
Frnax as & functiou of m in agreement with the numeri- 
cal solution of Fmax which is depicted in the inset of Fig. 
|. The results of F^nax and B{Fc) = minN C^ax / VN 
as a function of m show that the solutions of the two 
minimization problems are improved by increasing the 
number of phases m. 
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FIG. 3. minjv Cmax/V7^ vs. m as was obtained from the 
numerical solution of Eqs. and (^2|). The solid line is the 
least square fit 0.54/y^. 
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FIG. 4. Results for min ivC*max as a function of vn for se- 
quences of length = 32. The filled circles stand for the 
simulated annealing results. The deviation from the analyt- 
ical results (solid-line) is attributed to finite size efi'ects and 
to the suboptimal solution obtained in our limited running 
times of the simulations. Inset: Analytical results for Fraax 
as a function of m which were obtained from the numerical 
solution of Eqs. (^) and (|l2|). The solid line is the least 
square fit, 3.7m. 



These results raise the question whether it is possi- 
ble to increase m such that minjv Cmax becomes 1. An 
asymptotic expansion of Eqs. (|l^) and (|l2|) reveals that 
for m = N, min^v Cmax = 1- Note that C'n^i is always 1 
and therefore min^r Cmax = 1 holds for the entire regime 
TO > iV. 

We used the simulated annealing method for sequences 
of lengths = 32 with different number of phases to, to 
find miuAT Cmax- The relatively small sequences size was 
chosen to enable an appropriate scan of the configura- 
tion space in reasonable computational time. The re- 
sults are exhibited in Fig. ^, and support the anticipated 

N/m behavior of uAtin Cmax- For to = 32, there is 
a deviation of minjv Cmax from the analytical prediction 
min N Cmax — 1 , probably since the simulated annealing 
method yields only suboptimal solutions. 

Finally, we would like to examine the simulated an- 
nealing method in light of the results of this study. Sim- 
ulations show that for the same running times, the sim- 
ulated annealing method yields Cmax which is closer to 
its minimal value than X]^=i C^- This can be explained 
by the larger degeneracy of Cmax compared with that of 
S^=i Cx- Moreover, it turns out that in order to mini- 
mize Cmax , it is preferable to start the searching process 
with the minimization of the energy function, C^i 
and then replace it with Cmax [|l2|. In this way the sys- 
tem avoids the plateaus which characterize the landscape 
of Cmax in the configuration space. However, future re- 
search is necessary to find out how the results of this 
study can be applied to further improve the searching 
processes of low autocorrelated sequences. 
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